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ABSTRACT 

Context. Our position inside the Galaxy requires 3D-modelling to obtain the distribution of the Galactic magnetic field, cosmic-ray 
(CR) electrons and thermal electrons. 

Aims. Our intention is to find a Galactic 3D-model which agrees best with available radio observations. 

Methods. We constrain simulated all-sky maps in total intensity, linear polarization, and rotation measure (RM) by observations. For 
the simulated maps as a function of frequency we integrate in 15' wide cones the emission along the line of sight calculated from 
Galactic 3D-models. We test a number of large-scale magnetic field configurations and take the properties of the warm interstellar 
medium into account. 

Results. From a comparison of simulated and observed maps we are able to constrain the regular large-scale Galactic magnetic 
field in the disk and the halo of the Galaxy. The local regular field is 2 fiG and the average random field is about 3 fiG. The known 
local excess of synchrotron emission originating either from enhanced CR electrons or random magnetic fields is able to explain 
the observed high-latitude synchrotron emission. The thermal electron model (NE2001) in conjunction with a proper filling factor 
accounts for the observed optically thin thermal emission and low frequency absorption by optically thick emission. A coupling factor 
between thermal electrons and the random magnetic field component is proposed, which in addition to the small filling factor of 
thermal electrons increases small-scale RM fluctuations and thus accounts for the observed depolarization at 1.4 GHz. 
Conclusions. We conclude that an axisymmetric magnetic disk field configuration with reversals inside the solar circle fits available 
observations best. Out of the plane a strong toroidal magnetic field with different signs above and below the plane is needed to account 
for the observed high-latitude RMs. The large field strength is a consequence of the small thermal electron scale height of 1 kpc, which 
also limits the CR electron extent up to a height of 1 kpc not to contradict with the observed synchrotron emission out of the plane. Our 
preferred 3D-model fits the observed Galactic total intensity and polarized emission better than other models over a wide frequency 
range and also agrees with the observed RM from extragalactic sources. 

Key words, polarization - Radiation mechanisms: non-thermal - Radiation mechanisms: thermal - ISM: magnetic field - ISM: 
structure - Galaxy: structure 
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^-^ 1. Introduction maps. Polarization maps are influenced by Faraday rotation ef- 
fects along the line of sight and thus require to take the distri- 

Synchrotron radiation is the major component of the Galactic ra- ^^^^^^ ^f ^^^^^^^ electrons in the Galaxy into account as well, 

dio continuum emission from centimetre to longer wavelengths, important information on the Galactic magnetic field strength 

The diff^use Galactic synchrotron emission originates from rel- ^„^ direction comes from new rotation measure (RM) data. All 

H ativistic electrons gyrating in the magnetic field. An addi- ^^^^^ information must be commonly considered to constrain 

tional smaller contribution is observed from supernova remnants q^j^^^j^ 3D-models in order to find the most realistic one. This 

(SNRs) and pulsar wind nebulae located in the Galactic plane, p^p^^ p^g^^^j^ ^ gg^^^^l ^^^^^^^ ;„ jj^^j direction. 

where cosmic ray (CR) electrons are accelerated to high ener- ^^ . ,t^ , , ^ , • ^ , , ^t^ , 

e u » ■• ■ ■ » ■• 11 u- ui 1 • J u » Previous 3D-models tor the magnetic field, CR electrons and 

gies. Synchrotron emission is intrinsically highly polarized, but , , , , • , • T- • , ,, ,- 

c . ,. ,.• ff . • fu „ f- A ■ f f u ™ ^ „, thermal electrons were derived individually from selected data 

Faraday rotation errects in the magnetized interstellar medium „, . , , , ,,.,,, , . . 

, ;u 1- f • u» J 1 u • ff » J sets. Their general relevance can be established by combining 

along the line or sight and also beam averaging ettects reduce , , f , . , . . , \ . ° 

., u J \ 1 • »• » u 11 1 T» these models to explain other, in particular new observational 

theobservedpercentagepolarizationtomuchsmallervalues.lt , ^ . ,-„,.,.„., •, ,,,,^ 

is crucial to understand the regularity, strength and spatial struc- "^^'^^ !°' .'"^^^"'^^ ^'Tr^^ ^- "^f ^! '^°"''^^'' ^^' ^^^^ ^''" 

. f., » -c ij • /-^ 1 u- u »-ii J J- tragalactic sources (EGSs), thermal and non-thermal intensities 

ture or the magnetic field in our Galaxy, which are still under dis- ",. , , . , . . , ,. , 

uu u tu f: t J * t- e/^ i *• i • »• and linearly polarized emission together, which proves to give 
cussion although the first detection ot Galactic polarization was . ^ ^ . , ■ , ^%^ , , „ , , • 

, , ■ j „T- 1 , • — p — ~ni in£o L h\r * u -.ti l Stringent Constraints to Combined 3D-models. We note that sim- 

several decades ago (IWielebinski et al.1119621: IWesterho ut et al.l , ° , . , , ^ ^ , , • 

'T7JZ^^ /-v V- ,i.u/-'i *j-i »u ^oci J- mated maps are also required as a robust reference tor the techni- 

1962). Our position inside the Galactic disk at about 8.5 kpc dis- , ■ . ^^ , . ,., , r,,, . 



. e \u 1^ 1 ^- » • J 11- e ^u cal design of future telescope projects like the SKA or as inputs 

tance trom the Galactic centre requires modelling ot the mag- r , , ■ ■,.,■,, r , r^r »xt^t^ • ■ 

.. c ij J r^r> 1 » ^ ■ I * u 3 t, t_ tor data reduction pipelines like that tor the PLANCK mission, 

netic field and CR electrons to interpret observed synchrotron ^ ^ 

The paper is organized as follows: we review the status of 
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elling are discussed in Sect. 7. Conclusions are presented in Sect. 



2. Status of Galactic 3D-models 

Many attempts have been made to recover the structure of the 
Galactic synchrotron emission. Beuermann et al. ( 1985) decon- 
volve d the 408 MHz all-sky total intensity map dHaslam et alj 
11982 *) and modelled the emission as a thick disk with a kpc scale 
height contributing about 90% of the total emission plus an em- 
bedded thin disk with a scale height of abo ut 200 pc, which is 
consistent with the results bv Phillipps et al.l (1198 1 a b) using the 
same data. They adapt a spiral pattern for the emissivity and as- 
sume comparable regular an d random magnetic field strength. 
[Fleishman & TokarevI (Il995l) present observational evidence for 
an excess of synchrotron emissivity within a few hundred par- 
sec of the solar system based on low-frequency absorption data 
towards a number of large HII regions . A rec ent study of local 
Faraday screens by IWoUeben & ReichI (12004*) seems to confirm 
this result. Since local emission is also visible at high Galactic 
latitudes the thick disk assumption might be re-discussed. 

The observed spectrum of Galactic synchrotron emission 
shows complex spatial variations due to processes such as en- 
ergy losses, diffusion of CR electrons or variations of the mag- 
netic field strength, which makes a modelling of the spectral 
behavio ur rather difficult. For the northern sky iReich & ReichI 
(Il988allbl) found spectral variations {Ty oc \fi^ with Tv the bright- 
ness temperature at frequency v and p the temperature spectral 

index) in the range /? 3.1 to yS 2.5 between 408 MHz 

and 1420 MHz. A t lower frequencies the spectrum is flatter 
iRoger et al] d 19991) find spectral variations between /? ~ -2.55 
andyS 2.4 between 22 MHz and 408 MHz. At higher frequen- 
cies the synchrotron spectrum steepens, but the a mount is un- 
clear s o far. For the 1420 MHz and 22.8 GHz survevs Rei ch et all 
(|2004|) calculate steeper spectra towards higher latitudes than 
belo\ y 1420 MHz. Howev er, the 22.8 GHz (K-band) WMAP 
data dHinshaw et alj|2007l) may contain an un known fraction of 
anomalous dust emission JDavies et al.ll2006l) and thus the syn- 
chrotron spectrum might be steeper than the total emission spec- 
trum. For any modelling it implies that the assumption of a con- 
stant spectral index is a simplification. 

Synchrotron emission is generated by CR electrons in the 
Galactic magnetic field, which are considered as independent. 
We are unfortunately located in the disk of our Galaxy, which 
makes the observation of the large-scale field difficult. From ob- 
servations of other spiral gala xies it is known t hat the field basi- 
cally follows a spiral pattern dBecket al.l ll996V There are many 
methods to obtain inform ation about the Galactic ma gnetic field 
observationally (see e.g. iHan & Wielebinskil d2002l) ). A pow- 
erful method is the analysis of pulsar RMs, which indicates 
from recent data that the fi eld direction i s opposite in the arm 
and the inter-arm regions ( Han et al1 l2006). Howev er, RM data 
from EGSs can be inte rpreted with less reversa ls CBrown et alj 
l2003btlHaverkornetal. 2006; B rown et alj2007h . The local reg- 
ular field strength is about 2 //G and increases to abo ut 4 nG 
at abo ut 3 kpc from the Galactic centre according to H an et alj 
d2006l) . where no coupling between the ma gnetic field s t rength 
and the thermal electron density is assumed. IStrong et al.l d2000l) 
estimated the total local field to be about 6 juG. The Ga lactic ran- 
dom field component was estimated to be about 6 jiG dHan et alj 
[2006). This is about a factor of 2 to 3 larger than the regular field 
component. This ratio is larger than it was assumed in earlier 
models. 



The distribution of Galactic CR electrons has still to be 
established. The emission at radio frequencies originates from 
GeV electrons, where a power law distribution is a good approx- 
imation. SNRs are the acceleration sites and thus the sources 
of CR electrons. However, the propagation of CRs in our 
Galaxy is very complex and modelling elTorts like GALPROP 
dStrong & Moskal enkoll 19981) are required. 

The properties for the magnetized thermal interstellar 
medium are crucial for the interpretation of polarization observa- 
tions, because Faraday effects strongly modulate the polarization 
emission along the line of sight and make modelling of Galactic 
polarization quite challenging. At low frequencies Faraday ef- 
fects mi ght be strong e nough to entirely smear out distant polar- 
ization dSokoloff et alj |1998). For the Galactic large-scale ther- 
mal component the NE2001 model based on pulsar dispersion 
measures (Cordes & Lazio 2002) represents the best knowledge. 
However, additional Faraday effects come from HII regions and 
also from Faraday screens hosting strong regular fields, but a 
low thermal electr on density dWoUeben & Reichl2004tlSun et alJ 
120071 : lReichl2006h . and might influence the observations by a yet 
unknown amount. 



3. The method applied 

The observables are RM data, total and polarized intensity maps 
at different frequencies. Polarization angles and percentage po- 
larization can be used also for constraints. The goal is to ob- 
tain the thermal electron, CR electron and magnetic field 3D- 
distribution in the Galaxy which agrees best with the observa- 
tions. 

RMs depend on the regular magnetic field component and 
the thermal electron density along the line of sight. The perpen- 
dicular magnetic field component and the CR electron density 
determine the total intensity synchrotron emission. Observed po- 
larized intensity depends on the regular fraction of the magnetic 
field, but gets modified as a function of frequency by the line of 
sight magnetic field components and the thermal electron den- 
sity. To achieve a 'best model' an optimized quantity like x^ 
should be used to adjust the model parameter This is beyond 
the scope of the present attempt and also clearly needs more ob- 
servations than actually available for a well established model. 

We tackle the 3D-modelling problem by finding the model 
parameters in four steps as shown in Fig. [1] The modelling is 
based on the thermal electron density distribu tion as it is de- 
scribed by the widely accepted NE2001 model ( Cordes & Lazid 
2002). After considering a proper filling factor for the thermal 
electron distribution of the NE2001 model we are able to repro- 
duce the observed optically thin free-free emission and the effect 
of absorption from the optically thick gas. This indicates that the 
NE2001 model is valid in general, although a larger scale height 
is favoured by our modelling results as discussed later. From RM 
data of EGSs we are able to obtain the regular magnetic field 
direction and its strength separately for the Galactic plane and 
towards higher latitudes by keeping the thermal electron model 
parameters from NE2001 fixed. From the observed 408 MHz to- 
tal intensity and the 22.8 GHz polarized intensity we constrain 
the random magnetic field component and the large-scale distri- 
bution of CR electrons using the regular magnetic field compo- 
nent obtained in the step before. If the random field component 
is increased or decreased the CR electron density must be re- 
duced or enlarged accordingly to fit the total intensity emission. 
Both components must be matched to fit the observed polarized 
intensity. 
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Fig. 1. Overview on the different steps applied to find the Galactic 3D-emission components and their properties. Here f-f stands 
for free-free emission, ne for thermal electons, B (b) vectors for the regular (random) magnetic field component, and I (PI) for total 
(polarized) intensity. CGPS and SGPS mean the Canadian and the Southern Galactic Plane Survey, respectively. The blue color is 
for observational data, the pink for the model and the yellow for NE2001. 



The local excess of synchrotron emission is considered in 
our model, which is able to account for the observed high- 
latitude emission as well. At that point our model is able to rep- 
resent the total intensity synchrotron emission and the intrinsic 
polarization properties as observed at high frequencies and is 
in general agreement with the observed extragalactic RM dis- 
tribution. In the last step, the observed depolarization proper- 
ties at 1 .4 GHz are accounted for by introducing a coupling fac- 
tor between the thermal electron density and the magnetic field 
strength on small-scales in addition to the thermal filling factor 

4. Available observations constraining the 
3D-models 

To constrain the simulated maps from Galactic 3D-models we 
consider various observed all-sky data sets. These maps contain 
the integrated line of sight information of large-scale Galactic 
emission, but also show all kinds of local foreground structures 
and distinct source complexes, which are not the subject of mod- 
elling in our case. We note that it is not always clear how to dis- 
tinguish and to correctly separate the observed structures into lo- 
cal and global emission components. Unresolved strong sources 
can be easily identified and are not represented by the models. 
Large-scale local features like the giant radio loops with diame- 



ters up to 100° clearly show up as distinct feature in the total in- 
tensity and the polarization maps and are also not included in the 
model, although they are difficult to separate from the observed 
maps. The large-scale highly polarized "Fan region" is difficult 
to trace in total intensities, but certainly contributes to the ob- 
served emission. Distinct SNRs and large HII region complexes 
are concentrated in the Galactic plane and are excluded from 
the modelling, but their position is known from high-resolution 
surveys along the Galactic plane. Faraday screens of consider- 
able size exist in the interstellar medium, which are weak emit- 
ters, but reflect their existence by their large RMs by rotating 
the polarization angles of the background emission. Only a very 
densely sampled RM grid of sources or multi-frequency obser- 
vations of diffuse Galactic polarization can reveal such features. 
These data are not available, but their importance is widely ac- 
cepted. "Key s cience projects" for L OFAR and the future SKA 
were defined (iGaensler et al.l 120041) . which exactly have these 
aims. Faraday screens might influence our modelling results to 
an unknown extent. All the uncertainties in the recognition of 
local or distinct features limits any Galactic modelling attempt 
to the reconstruction of the large-scale Galactic structures. 
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Fig. 2. The WMAP 22.8 GHz all-sky polarized intensity map 
(upper panel) and the 1 .4 GHz all-sky polarized intensity map 
(lower panel). The polarized intensities are shown greyscale 
coded from to 100 juK for 22.8 GHz and from to 570 mK 
for 1.4 GHz. 



4. 1 . All-sky total intensity maps 

All-sky total intensity surveys up to a few GHz, which are 
adjusted to an absolute zero-level, are essential to obtain the 
required large-scale knowledge of the Galactic synchrotron 
emission. Carrying out such surveys is rather time-consuming 
and therefore these surveys have relatively low angular reso- 
lution. Available all-sky surveys have been reviewed by Reich 
(l2003h . Among these the 40 8 MHz surve y (Haslam et al. 1983) 
and the 1420 MH z survey (lReichlll982l: iReich & Reich, 1981 
iReich et al.l200lh have the highest sensitivity and an angular res- 
olution down to 0°6 and thus are well suited for our purpose. 

Galactic radio emission at frequencies below a few GHz is 
dominated by synchrotron emission, in particularly out of the 
plane, relative to free-free emission. Thermal dust emission and 
probably spinning dust emission are becoming important for fre- 
quencies higher than 5 GHz or 10 GHz. The 408 MHz total in- 
tensity map is often used as a "synchrotron template" and in fact 
significant thermal emission is confined to the Galactic plane. At 
1420 MHz the thermal fraction in the plane increases t o about 
40% in the inner Galaxy (Reich & Reich 1988b). Paladi ni et all 
(l2005h estimated a thermal contribution with a maximum of 
about 68% for the region of 20° < Z < 30° and \b\ < 1?5, where 
discrete HII regions contribute about 9%. Note that these esti- 
mates are based on the separation of thermal and non-thermal 
components assuming appropriate spectral indices. However, the 
spectral index of synchrotron emission is uncertain and might 
vary with direction, although these estimates are largely consis- 



tent with the method to separate both compone nts based on a re- 
lation between infrared and thermal emission dBroadbent et al.l 
1989). It is certainly worth to refine the separation of thermal 
and non-thermal emission, but this is beyond the scope of the 
pre sent paper. 

[Dickinson et al.l (l2003h have constructed an extinction cor- 
rected free-free emission template for about 95% of the sky 
based on Ha surveys of the northern sky ( Ha ffner et al.l 120031) 
and of the southern sky ( Gaustad et al. 2001) by using dust emis- 
sion data ( 'Schlege l et al.il99 8) to calculate the extinction correc- 
tion. iHins haw et alj (l2007h used the Ha template by Fink beinen 
(12003 ') as a prior estimate in the maximum entropy method anal- 
ysis of the 3-year WMAP data to construct the free-free emis- 
sion template. This template might contain some leakage from 
other components such as synchrotron emission and anoma- 
lous dust emission. However, actually it probably represents 
the best estimate for low latitude emission. W e will refer to 
the free-free emission map at 22.8 GHz from iHinshaw et al.l 
( 20 0j) for comparison with the modelled thermal emission. The 
Dickinson et al.l ( l2003h template cannot be used for the Galactic 
plane area because of the large and uncertain extinction there. 

All radio maps show a clear concentration of emission along 
the Galactic plane. As mentioned above, local structures out of 
the plane are rather dominant features. The most striking pro- 
trusion emerging at / ~ 30° out of the plane is called the North 
Polar Spur (NPS), which can be traced up to very high latitudes. 
This spur is a part of the giant Loop I, which is believed to be 
a local old supern ova remnant at a distance of about 200 pc 
(lBerkhuiisenlll973l) . A recent new model by IWolleben I (|2007|) 
tries to explain the NPS in total intensity and linear polarization 
by local synchrotron emitting shells with a swept-up magnetic 
fie ld. A number of weaker lo ops of similar size exist as discussed 
bv lBerkhuiisen et alj ( Il971h . 



4.2. All-sky polarization maps 

The status o f Gala c tic po larization surveys has been recently 
reviewed by ' ReichI (12006'). For the interpretation of polariza- 
tion surveys the correct absolute zero -level is of particul ar im- 
portance as it was demonstrated by iReich et all (|2004|) . The 
only ground based all-sky polarization map exists at 1.4 GHz 
(Reich et al., in prep.) combining the no rthern sky survey car - 
ried out with the DRAO 26-m telescope flWolleben et al.ll2006l) . 
whi ch is tied to th e absolute level of the 1 .4 GHz Leiden sur- 
vey dBrouw & Spo elstra 1976), and the Villa Elisa southern sky 
survey ("Testori etal.ii20 04, Testori et al., in prep.), which has a 
large overlap with the DRAO survey. At this frequency the depo- 
larization effects are already very strong at latitudes below about 
30°, especially towards the Galactic plane as can be seen from 
Fig.|2] Thus the distance up to which polarized emission can be 
observed at 1 .4 GHz with an angular resolution of 36' is small 
compared to the line of sight through the Galaxy. O n the other 
hand the WMAP 22.8 GHz polarized intensity map (Pag e et al] 
2007) is almost unaff'ected by Faraday depolarization effects and 
therefore should represent the intrinsic synchrotron properties 
of the Galaxy, since dust polarization is believed to be entirely 
unimportant at that frequency (Page et al. 2007). 

As it can be seen from the 22.8 GHz map there are strong im- 
prints of Loop I, which shows strong polarization along the NPS 
and also below the Galactic plane. Loop I is very local as men- 
tioned in the previous section and therefore can be also traced 
for large sections at 1.4 GHz. Another very bright polarization 
feature called the "Fan region" is centered at Galactic longitude 
of about 140° and was recognized as a local feature at a distance 
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of ab out 500 pc dWilkinson & Smithll 19741: iBrouw & Spoelstral 
11976 *). However, a new analysis based on the DRAO 1.4 GHz 
survey suggests that it migh t extend from about 500 pc to about 
2 kpc JWoUeben et al.ll2()()6l) . We take both structures as distinct 
local features, as they are very bright at 1.4 GHz contrary to 
the polarized emission from the Galactic plane. The modelling 
of these distinct features is beyond the scope of our simulation. 
They will not influence our modelling of the large-scale Galactic 
structures. 

Although external Faraday depolarization is very small at 
22.8 GHz we find a low percentage polarization at 22.8 GHz 
in gen eral. This is in agreement with the results of iKogut et alj 
( l2007h . We divided the polarized intensity map by the syn- 
chrotron emission template provided by the WMAP team. The 
polarization percentage in the Galactic plane is typically below 
2%. At higher latitudes typical values between 12% and 30% are 
derived, which are rather similar to those found at 1.4 GHz. Of 
course, the 22.8 GHz percentage polarizations may be system- 
atically underestimated in case of the existence of a significant 
spinning dust component in the total intensity map. Despite of 
this uncertainty the intrinsic synchrotron polarization will stay 
significantly below the maximum polarization value for a per- 
fect regular field. This might indicate that the random magnetic 
field exceeds the regular magnetic field component. Otherwise 
the superposition of a number of regular magnetic field compo- 
nents along the line of sight with different orientations may have 
the same effect. 

Comparing the 22.8 GHz and the 1.4 GHz polarized emis- 
sion we see strong depolarization eff'ects along the Galactic 
plane except for the "Fan region". This depolarization is at- 
tributed to Faraday effects from the warm ionized medium. 
However, it is not immediately clear, whether the diffuse 
medium or sources like HII regions or Faraday screens play the 
main role in this process. We will show below that the depo- 
larization calculated from the diffuse thermal component as pro- 
vided by the NE2001 model is insufficient to explain the 1 .4 GHz 
polarization observations. 

The polarized intensity can be used to constrain the strength 
of the regular magnetic field. As can be seen from the 22.8 GHz 
polarization map (Fig.|2l) the main characteristics are (1) strong 
polarization towards the Galactic centre; (2) the distribution of 
polarized emission along Galactic longitude is asymmetric with 
respect to the Galactic centre. Depolarization can be neglected 
at 22.8 GHz and therefore the polarization angles could also be 
used to check the orientation of the magnetic field. However, 
our simulations show that any spiral magnetic field configuration 
gives rather similar polarization angle patterns. Therefore we do 
not use polarization angles to decide on the fi eld configuration 
as it was previously done bv lPage et alj ( l2007h . As shown below 
this model is inconsistent with other observations. 

4.3. Rotation measures of extragalactic sources 

When a linearly polarized electromagnetic wave passes a mag- 
netized thermal medium its polarization angle experiences rota- 
tion, called Faraday rotation. The rotation is proportional to the 
square of the wavelength and the coefficient is the RM, which 
can be written as 




VK' rad/B"^ 



RM^ 0.81 n,(l)Bn(l)dl 



(1) 



where the integral range is along the line of sight from the ob- 
server to the source at a distance D measured in pc. n^ per cm""* 



Fig. 3. Synthesized RM map constructed from published and 
new RM data of EGSs (Han et al., in prep.). The preliminary 
new RM data set from the Effelsberg survey was provided by 
JinLin Han. The RM map does not include Galactic plane RM 
data from the Canadian (CGPS) or the Southern Galactic Plane 
(SGPS) surveys. 



is the thermal electron density and B\\ in fiG is the magnetic field 
component along the line of sight in the intervening medium. 
RM is positive when the magnetic field points towards us. 

RMs of EGSs reflect the integrated effect of the magnetized 
medium throughout the Galaxy and can be used for a diagnosis 
of the thermal electron density and magnetic field strength along 
the line of sight. The RM contribution from outside the Galaxy 
should be very small. The effect of the EGS-intrinsic RM can 
only be accounted for by averaging a sufficient large number 
of RMs in a small area. Unfortunately RM data of EGSs are 
sparsely distributed in the sky despite recent observational ef- 
forts. This holds in particular for the southern sky, which makes 
the construction of an all- sky RM ternplate a difficult task. With 
the RM data compiled bv lHan et alj d 19991) and the preliminary 
RM data from the recent L-band Effelsberg RM survey of about 
1800 polarized NVSS sources (Han et al., in prep.) we synthe- 
sized the RM all-sky map show n in Fig. [3l following the same 
interpolation method as used by Uvaniker et alj (11998.) to create 
a smooth map from ungridded data. The observed data close to 
each pixel of the interpolated all-sky RM map are added with 
a weight calculated from the function exp(-(0/0)™), where 6 is 
the angular distance to the pixel. was set to 8° and m to 0.3. 
This R M m ap contains more data than us ed bv lDineen & ColesI 
(|2005[) and IJohnston-HoUitt et al.l (|2004|) . However, flie basic 
large-scale features seen in all RM maps are rather similar. 
Although there are less data in the southern sky, which causes 
edge effects as visible in Fig. [3] the data clearly indicate the gen- 
eral RM trend as discussed be l ow. We did not include the numer- 
ous RM data bv lBrown et al.l (l2003al 12007) from sources in the 
Galactic plane as the intention is to trace the high-latitude RM 
characteristics. We treat the Galactic plane RM characteristics 
separately. 

The RM map in general shows an antisymmetric distribu- 
tion below and above the Galactic plane an d also with respec t to 
Galactic longitude 0°, as already noted bv iHan et al.l (Il997h . In 
a recent RM study towards the Galactic centre region thi s RM 
pattern was also seen towards the inner Galaxy ( lLawll2007l) . The 
high-latitude RM values are about 20-30 rad m"^ on average. 

For the RM distribution along the Galactic plane we rely 
on the EGSs detected in high-resolution Galactic plane surveys, 
where sources contrast well against the Galactic diffuse emis- 
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Fig. 4. RM profiles along the Galac tic plane with RM da ta for the 
northern Galactic plane taken from lBrown e t al.' ('2003a") and RM 
data for the southern Galactic plane from Brown et al. (20071) . 
The average for each longitude bin is shown in the lower panel. 
The bars indicate the +lcr variance. 



sion. The Canadian Galactic plane survey (CGPS, iTavlor et al] 
l2003h covers parts of the northern sky. The south ern Galactic 
plane is covered by the SGPS dG aensler et all bOOl ). The 
RM data from the CGPS dBrown eT al. 2003^ and the SGPS 
dBrown et alj2007l ) are shown in Fig.|4] Note that the fluctuation 
of RMs is very large, typically around several hundred rad m"^, 
and averaging over large areas is needed to trace systematic RM 
changes. We sort the data into longitude bins and obtain their 
average to be compared with the simulations. F or the SGPS 
data th e binning procedure is the same as used by i Brown et alj 
(120071) . For the CGPS data, the RMs between I = 82° and / = 96° 
are divided into three bins containing 38 or 39 sources and the 
RMs elsewhere are put into 1 1 bins containing 37 or 38 sources. 
The binned data are also shown in Fig. |4] We note the follow- 
ing RM structural features: (1) the average RMs run smoothly 



from about -100 rad m"^ at I ^ 140° to about -200 rad m"^ at 
/ ~ 100°; (2) the sign of the RMs changes at / x 304°. We beheve 
these features to be real. The local minimum of RMs at Z » 310° 
seems also reliable. There are some more features, which seem 
less certain like RM values close to zero appearing at I » 320° 
and / ^ 275°. The RM fluctuations in the Galactic centre re- 
gion are extremely large (iRov et alj|2005l) . which seem to reflect 
strong variations of the properties of the interstellar medium lo- 
cated in the centre region. The available amount of RM data is 
insufficient to allow to constrain the large-scale properties of the 
interstellar medium and the magnetic field in the central region 
of the Galaxy. 

The large RM fluctuations observed might originate from 
the intrinsic RMs of the EGSs or the interstellar medium in the 
Galaxy. Fluctuations of the Galactic magnetic field have been 
estimated to be about 6 juG CHan et al.ii2004.) . which might cause 
the large scatter of RMs. However, we will show below that the 
dumpiness of thermal electrons and a likely correlation with the 
random magnetic field component causes large enough RM fluc- 
tuations for depolarization in accordance with the observations. 
We note that strong RM fluctuations are a handicap for the detec- 
tion of cosmological magnetic fields from fluc tuations of fore- 
ground corrected RMs as reviewed bv lWidrowl (120021) . 

Currently the RMs of pulsars are not considered for mod- 
elling. First of all the distances of pulsars are quite uncertain as 
many of them are obtained from the dispersion measure based on 



global electron density models, e.g. NE2001. Second, our sim- 
ulation results are all-sky maps, which are difficult to compare 
with the RMs of pulsars, which are distributed along the line 
of sight. Finally, despite of recent efforts, pulsar RM data are 
still sparsely distributed in the Galactic plane and it is difficult 
to account for the influence of the local interstellar medium. In 
future the pulsar RM data definitely will be used as important 
additional constraints in particular as they have no intrinsic RM 
contribution. 
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Fig. 5. The NE2001 thermal electron density model by 
ICordes & Laziol ( l2002l |2003|) projected on the Galactic plane. 
The position of the Sun as well as the designation of the spiral 
arms are indicated. 



5. Numerical realization of Galactic 3D-modelling 



IWaeUcensI (l2005h developed the HAMMURABI code aiming to 
simulate the all-sky synchrotron emission of the Galaxy. The 
code takes Faraday rotation effects into account and thus Faraday 
depolarization effects are naturally involved. In this work we im- 
plemented additional entries to account for optically thick free- 
free emission. Th e code output map s follow the HEALPIX pix- 
elization scheme dGorski et alj|2005l) . which divides the sky into 
pixels of equal areas. The total number of pixels for an all-sky 
map, which defines its angular resolution, is A^pix - 12A^gjpg, 



with A^siDE = 2*= and k = 0,1,2,3, 
pixel (A6) can be approximated as. 



A6I: 



3 3600' 

TT A^SIDE 



The angular size of a 



(2) 



In our simulation we use A^side = 256, which corresponds to 
an angular resolution of 15'. This selection is a compromise be- 
tween angular resolution and computing time. 

Along the line of sight the cone for each pixel is divided into 
volume elements in radial direction, where in this work each unit 
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has an equal radial length Ar = 20pc. The largest volume unit in 
such a cone is approximately 

Ar 

(3) 



''max — ^•'"max 



12-A^de' 



where r^jax is the longest possible line of sight with the maxi- 
mal value of 32 kpc. Due to the disc like shape of our galaxy 
the effective distance out to which we perform our computations 
varies and decreases sharply towards high latitudes and towards 
the anti-centre direction. Note that fluctuations of the magnetic 
fields on scales smaller than the volume unit cannot be resolved 
by the code. 

For the /-th volume element of a cone we compute the fol- 
lowing quantities, 

Pli = Cp,B|_'7'^^(i+'"/2^r 
RMi = 0.81n,,B,-,||Ar 

Ui = Plism(2i//i) (4) 

Qi = P/,cos(2tA,) 
EMj = n^,.Ar 

Ti = 8.235 X 10-2rri-35v-2i£M; 
, If = r,(l-exp(-T,)) 

Note that the calculated Stokes values Uj and Qi for each vol- 
ume element include the effect of foreground RM. The CR elec- 
trons are assumed to follow a power law with an energy spectral 
index p. The values C/ and Cpj depend on the spectral index 
p and are further detailed in iRybicki & Lightnianl (Il979l) . The 
synchrotron emission is related to the magnetic field component 
perpendicular to the line of sight B,, ^, while the RM depends on 
the magnetic field parallel to the line of sight B, y. The intrinsic 
polarization angle in each volume element is i^,; o, which is de- 
fined as the inclination angle of Bjj^ with respect to the north 
(in the frame of Galactic coordinates). Note that the magnetic 
field is the sum of the regular and random field components. The 
thermal electron density is n^,, and its temperature is T,. For the 
re lation betwe en opacity r, and emission measure, EM, we refer 
to lRohlfs& W ilson (2000). The calculation of EMi, r, and if 
are expansions of the original HAMMURABI code. 

Intensities and RMs for each pixel are obtained as the inte- 
gral of the contributions from all the volume units along the line 
of sight 

Q = E,e/ (5) 

u = ZiUi 

[RM = ZiRMi 

All intensities are calculated as a function of frequency. 
More details on the HAMMURABI code will be discussed in 
a forthcoming paper by WaeUcens, EnBlin, et al. (in prep.). 

6. Modelling 

For the description of the 3D-models we introduce cylindrical 
coordinates (R, (p, z), where R is the Galactocentric radius, is 



the azimuth angle starting from / - 180° and increasing in anti- 
clockwise direction, and z is the distance to the Galactic plane. 
To show the slices of the 3D models a Galactocentric Cartesian 
coordinate is also used, where z is the same as in the cylindrical 
system, and the x-y plane coincides to the Galactic plane with x 
pointing to / = 0° and y to I - 90°. Through all calculations the 
distance of the Sun to the Galactic centre is taken as 8.5 kpc. 

Below we present a detailed 3D modelling of the Galactic 
emission following the steps described in Sect. 3. 

6.1. Thermal electron density 

6.1 .1 . The model for diffuse ionized gas - NE2001 



According to iFerriera (1200 ih the warm ionized medium (WIM) 
has an electron density of 0.2 cm"^ to 0.5 cm"^, a volume filling 
factor of 5% to 14% and a temperature of about 8000 K. It con- 
sists of HII regions embedded in diffuse ionized gas (DIG). The 
average electron density of HII regions can be derived from its 
radio continuum integrated flux density (e.g. ISchraml & Mezgen 
Il969h . However, the distances of HII regions are required, which 
are for instance derived from observations of Doppler shifted re- 
combination lines. Kinematic distances are based on Galactic ro- 
tation curves and are not quite certain. The scale height for HII 
regions is a bout 50 pc in the inner Galaxy and increases towards 
larger radii (iPaladini et al. 2004). 

HII regions along the line of sight can influen ce RMs from 
sources in the background as shown bv lMitra et al.l ([2003.) . Since 
HII regions are tightly confined to the Galactic plane, the chance 
that a line of sight passes an HII region outside the plane is small. 
In the Galactic plane the typical effect of HII regions on RMs is 
difficult to estimate in general. Both the electron density and the 
magnetic field inside a HII region may be an order of magnitude 
larger compared to the diffuse gas. A t ypical RM c ontribution of 
several tens of rad m ^ was quoted bv lMitra etal.1 (120031) . These 
values are smaller compared to the average RM of EGs, which 
trace the entire magnetized DIG in the Galaxy along the line of 
sight. Therefore we consider only the effect of the DIG in our 
simulations. 

A good tracer of the DIG is the dispersion measure (DM) 
of pulsars. The DM is defined as the integral of the elec- 
tron density along the line of sight. Using pulsars with known 
distances as determined from e.g. parallax or Hl-absorption 
measurements an electron density model can be obtained by 
a fit of the observed DMs. In this wa y a two-disk mode l 
was proposed and devel oped in detail bv iCordes et al.l (0*99 lb; 



iTavlor & Cord es ( 1993 ): iGomez etal.1 (|2001|) . ICordes & Laziol 
(l2002h and ICordes & Laziol (1200 3*) used more pulsars and im- 
proved the model by involving local irregularities and in par- 
ticular added a spiral arm structure. Thus the NE2001 model has 
three components: a thick disk, a thin disk and spiral arms as dis- 
played in Fig.|5] Locally the mid-plane electron density is about 
0.03 cm"^ and the exponential scale height is about 0.97 kpc. 



6.1 .2. Tine filling factor of the DIG - from free-free emission 
at 22.8 GHz 

The properties of the DIG can also be traced by its free-free 
emission. This is done either by measuring Ho- intensities or ra- 
dio continuum emission. However, the diffuse Galactic thermal 
and non-thermal radio emission components mix and are diffi- 
cult to separate. Ho- observations need an extinction correction t o 
be a reliable free-free emission template dDickinson et al.l2003h . 
which is, however, only possible out of the plane. Because the 
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Fig. 6. Longitude and latitude profiles from the simulated and the free-free emission maps at 22.8 GHz. The dotted lines indicate 
the data from the 22.8 GHz template. The solid and dashed lines are from the thermal emission simulations based on the NE2001 
model without (dashed) and with (solid) a filling factor correction according to lBerkhuiisen et al., (2006). 



distances of the emission components along the line of sight are 
unknown, the electron density distribution cannot be obtained di- 
rectly. However, the free-free emission template can be used as 
a check to prove that the integrated properties of the DIG model 
are appropriate. 

Before we calculate the all-sky free-free emission map from 
the NE2001 model we assess the contributions from HII regions. 
The thermal emission from HII re gions is confined to a thin layer 
{\b\ < r e.g. Paladiniet all (12004)) along the Galactic plane with 
a small scale height. We do not include individual HII regions in 
our simulations of the free-free emission, which we believe is an 
acceptable simplification. 

The electron temperature is required to calculate the free- 
free emission intensity as seen from Eq. (|4|. The existence of an 
electron temperature gradient in the Gala xy is known and was 
recently confirmed with some precision bv lOuireza et alj (l2006l) 
from recombination line measurements as Ti(R) = (5780±350)-i- 
(287 + 46)R, where R is in units of kpc. This corresponds to 
an electron temperature of about 800 K near the solar s ystem, 
which is larger than that derived by iPaladini et al.l ( l2004l) . Note 
that these temperatures are obtained from observations of HII 
regions, which are k nown to be l ower than the DIG by typically 
about 1000-2000 K dMadsen et al. 2006) . However, these tem- 
perature differences have just a small influence on the amount 
of free-free emission (see Eq. (|4|). For simplicity we use the 
temperature gradient established from the measurements of HII 
reg ions. 

[Reynolds et al.l (Il999h reported a tempe rature increase with 
increa sing z- This increase was fitted by [Peterson & Webben 
(120021) as Ti(z) = 7000 - 526z + 1770z^ where z is in units of 
kpc. Combining both temperature gradients we obtain !',(/?, z) - 
(5780±350)+(287±46)/?-526z+ 1770^2 being used in our calcu- 
lation. Although the temperature is high at large z the free-free 



emission remains small due to the low electron density there. 
Our results are not sensitive to the temperature gradient in z- 
We note that the ele ctron gradien t in z direction is still debated. 
Recent results from iDavies et alj ( 1200 6) based on WMAP and 
ancillary data indicate rather low electron temperatures at high 
latitudes. 

When we calculate the free-free emission from the NE2001 
model as shown in Fig. |5] by using Eq. (HJi, we obtain much 
smaller intensities than from the optically th in free-free emis- 
sion template derived bv lHinshaw et alj (l2007h based on WMAP 
data. This can be seen from Fig. |6] where the longitude pro- 
files (average for \b\ < 5°) and the latitude profiles (aver- 
age for the longitude ranges indicated in the plots) are shown. 
However, the differences vanish when the clumpy nature of the 
DIG is taken properly into account, which means to introduce 
an appropriate filling f actor of the thermal electrons. Recently 
Berkhuii sen et al.l 12006) studied the Galactic filling factor /^ of 
the DIG and its dependence on z in detail. The filling factor /^ is 
defi ned as fe =< rie >^ I < nj >, where the <> means the aver- 
age. |Be^^isenitaD (iOO^ derived /,(z) = 0.07 exp(|z|/0.5), 
where z is in units of kpc, and /« = 0.32 for z larger than 
0.75 kpc. EM calculated from the NE2001 model is based on 
< lie > and needs to be corrected by the factor life, which 
means that EM increases by a factor 11 fe- After employing the 
filling factor correction the comparison between the free-free 
emission template from WMAP observations and the corrected 
NE2001 model is shown in Figs. |6] and |2l The simulations now 
resemble the WMAP template quite well. From the comparison 
shown in Fig. |6]it is also obvious that the thermal high-latitude 
emission from the simulations is slightly higher than that from 
the template, which might indicate a further increase of the fill- 
ing factor towards larger z beyond the limit of 0.32 assumed by 
iBerkhuiisen et al.1 (l2006h . 
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Fig. 7 . The WMAP free-free emission template (iHinshaw et alJ 
120071) is shown in the upper panel and the NE2001 based simu- 
lated maps after taking a filling factor into account is shown in 
the lower panel. Both maps are shown at an angular resolution 
of 2°. 
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6.1 .3. The filling factor of the DIG - from low-frequency 
absorption 



Fig. 8. The average spectral index along the Galactic plane 
{\b\ < 3°) calculated from two simulated maps at 408 MHz and 
22 MHz including thermal absorption from the NE2001 model 
with (solid line) and without a filling factor correction (dashed 
line). For comparison we show the observed spectral indices 
(1^1 < 3°) calculate d from the 22 MHz nort hern sky survey map 
(iRogeretalJI 19991) and the 408 MHz map (lHaslametalJll982h 
at an angular resolution of 2° (circles). 



6.1.4. Remarks 

The NE2001 model represents the average large-scale electron 
density distribution, including some sources and large com- 
plexes like the Gum nebula. However, some similarly prominent 
features like the Cygnus complex at about I ~ 90° are missing 
as it is clearly seen in the free-free emission template and the 
low-frequency absorption spectra. This seems to be related to 
the much larger number of pulsars in the southern sky, which 
allows a modelling of such features there. After considering ap- 
propriate filling factors the NE2001 model is able to reproduce 
the WMAP free-free emission template at 22.8 GHz and also the 
low-frequency absorption quite well. 

In summary th e NE2001 model combin ed with the filling 
factor obtained by iBerkhuiisen et al.l (l2006l) are proven to be 
sufficient to model the low latitude DIG. We will discuss later 
that there are observational problems suggesting that at high- 
latitudes the thermal electrons density is underestimated. 



At low frequencies the optical thickness of free-free emission in- 
creases (Eq.m which implies that the synchrotron emission from 
behind is attenuated by a factor exp(-T). This can be seen by an 
apparent increase of the spectral index by t/ ln(v/,/v;), where v/, 
and v/ are the high and the low frequency, respectively. Since the 
optical thickness is proportional to EM the filling factor correc- 
tion is essential. We simulate total intensity maps at 408 MHz 
and 22 MH z by using an obser ved spectral index of fi - -2.47 
obtained bv lRoger et al] (IT999I) from the 408 MHz and 22 MHz 
surveys for regions outside the Galactic plane, which closely 
represents the spectral index of the Galactic low frequency syn- 
chrotron emission. The spectral index map is calculated from the 
simulated maps at the two frequencies and the average profile 
along the Galactic plane is compared with that observed (Fig.|8]i. 
We find that the spectral index increase from / = 60° towards the 
Galactic centre is well reproduced by our simulation. Note that 
the peaks are from individual source complexes, which are not 
included in our model. For comparison we show also the spectral 
index profile from the NE2001 model without the filling factor 
correction, which does not show the spectral increase towards 
the inner Galaxy caused by thermal absorption. 



6.2. Regular magnetic field properties 

The Galactic magnetic field is commonly described as a vector 
with a random component b and a regular component B. It can 
be written as, Btot = B + b. Following the procedure described in 
Sect. 3, the random field component will be discussed together 
with the CR electron distribution. The regular field properties are 
constrained by RMs. Below we use B = |B| to denote the strength 
of the magnetic field vector We treat the regular field component 
as the combination of a disk field B^ and a halo field component 
B^. The former is constrained by RMs near the Galactic plane 
(\b\ < 5°) and the latter by an all-sky RM map from EGS data. 



6.2.1. The disk field 

It is generally accepted that our Galaxy has an organized large- 
scale di sk field similar to other n earby Galaxies (s ee reviews 
by e^. iHan & Wielebinskil 120021) . As described by iBeck et d] 
(199a) the field more or less follows a spiral pattern that can 
be basically classified into an axi-symmetric spiral (ASS) with 
no dependence on the azimuthal angle or a bi-symmetric spiral 
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(BSS), which has a symmetry of tt. Relative to the Galactic plane 
the field above and below can be either symmetric (even parity 
or quadrupole) or asymmetric (odd parity or dipole). 

Magnetic fields are often discussed in the context of prop- 
agation of high-energy cosmic rays. The ASS field configura- 
tion has bee n use d by Stanev (199 7i) and the BSS field con- 
figuration by 'Stanev ( 19971) . iTinvakov & TkachevI (l2002h and 
iProuza & Smida (2003i) for modelling. A logarithmic spiral arm 
field was proposed bv lPage et al.l ( l2007h to interpret the WMAP 
observations. We collected a number of published models and 
checked them for their ability to simulate the observed RMs in 
the Galactic plane. The results are shown in Fig. |9] Obviously 
none of the models is able to reproduce the observed systematic 
features. This motivates us to revise the parameters previously 
used. 

The general form of a coplanar and constant pitch angle spi- 
ral magnetic field configuration can be written in the cylindrical 
coordinate system as 






DiiR,4>,z)D2iR,(p,z) sin p 
B^ = -Di(R,cf>,z)D2(R,cl>,z) cos p 
B^ ^ 



(6) 



where Di(R,(p,z} constrains the spatial variation of the field 
strength and D2(R,(p,z) introduces reversals or assymmetries. 
The definition of the pitch angle p is the same as previously used 
bylHan & Oiao ( 1994). We use the function Di(i?, (p, z) always in 
the following way: 



Dx{R,z) 



B„exp|-^-Mu>7?. 



(7) 



Br 



R<Rc 



However, for different models the variables Rq, Rc and zo 
may change. Th is kind of variation is the same as used by 
iHan et all (12006') to describe the magnetic field from pulsar RM 
observations. Too few RM data were observed so far in the first 
quadrant, which does not allow a detailed modelling of the field 
within the inner Galaxy. Therefore we keep the field strength 
constant there. 

The regular z-component of the Galactic magnetic field is 
ve ry small. Its strengt h is about 0.2 fiG to 0.3 fiG as estimated 
bv lHan & Qiao! (|1994|) . We assume B? to be zero in our models. 

Below we list three different magnetic field models: (1) the 
ASS model plus reversals in rings (ASS-I-RING); (2) the ASS 
model plus reversals following arms (ASS+ARM); and (3) the 
BSS model. Following the procedure described in Sect. 3 we 
use the NE2001 model for the thermal electrons to adjust the 
parameters in the three models by a trial-and-error approach in 
order to achieve the best fit for the observed RMs in the plane. 
The results are summarized below: 

ASS+RING — In this model, the parameters in Di{R,z) are 
Rq = 10 kpc, zq - I kpc, Rc - 5 kpc, Bq = 2 jjG and Be -2 jjG. 
The field reversals are specified by D2{R) being symmetric in z 
as. 



D2(R) 



+ 1 R> 7.5 kpc 

-1 6 kpc < 7? < 7.5 kpc 

-1-1 5 kpc < R <6 kpc 

-1 R<5 kpc 



(8) 



Here -i-l means a clockwise direction as seen from the north 
pole. The pitch angle is taken as -12°, which is the average of 
pitch angles used for the spiral arms in the NE2001 model. 



It can be seen from Fig. [10] that the simulated RMs in the 
Galactic plane fit the data of the CGPS and the SGPS well. We 
also used the RM data from Rov et al. (2005) near the centre and 
elsewhere from the compilation of iHan et al.l (ll999J) . Although 
the number of measurements is small they basically support the 
results from the simulations. 

ASS+ARM — In this model the parameters for Di{R,z) are 
Rq - 8.5 kpc, Zo = 1 kpc, Rc - 5.3 kpc, Bp - 2 ^G and 
Be - 2 fiG. We use similar logarithmic arms as lCordes & Lazid 
(l2003l) (see Fig. [S]), but with a pitch angle of -12°. The pitch 
angle of the magnetic field is the same as that of the arms. The 
reversals are placed between the inner edge of the Sagittarius- 
Carina arm and the inner edge of the S cutum-Crux arm. T his re- 
versal is consistent with that found bv lBrown et alj (|2007|) based 
on the SGPS RM data. However, this should not be taken as a fi- 
nal conclusion because the amount of RM data is small and also 
the position and the extent of the spiral arms is uncertain. As can 
be seen from Fig. [TO] this model fits the observations as well as 
the ASS+RING model. 

BSS — In this model the parameters for Di{R,z) are Rq = 
6 kpc, Zo = 1 kpc. Re = 3 kpc, Bq = 2 /iG and Be - 2 juG. 
D2{R, (p) is written as 



D2{R, 0) = - cos 



R 

Rh 



(9) 



where/? - 1/ tan/?. The parameters areRh - 9 kpc and/? = -10° 
for R > 6 kpc and otherwise Rj, - 6 kpc and p - -15°. The 
model is able to fit the RM data from the SGPS and the CGPS 
(Fig. [Toll, but is inconsistent with other data as discussed below. 

6.2.2. The halo field 

The synthesized all-sky RM map (Fig. [3]l shows an asymme- 
try in longitude and latitude relative to the Galactic plane and 
the Galactic centre, respectively, indicating an asymmetric halo 
magnetic field configuration. 

Such a RM configuration could in principle also originate 
from large local shells. However, at least two shells with op- 
posite magnetic fi eld orientat i ons lo cated below and above the 
plane are needed. [Han et all d 19991) showed that the RMs of 
pulsars with distances larger than 1 kpc increase with distance, 
which makes an explanation by local shells unlikely and strongly 
supports the large-scale field scenario. 

We follow iProuza & Smidal ( l2003l) for a description of the 
regular double-torus field component for one half of the Galaxy 
and the reversed direction in the other half: 



B^(R,z) = B^- 



R 



1 + 



\z\ - z, 



H^-^'"P 



R-R^,\ 



(10) 



The parameters are z^ = 1 .5 kpc, zf = 0.2 kpc for \z\ < Zq and 
otherwise zf = 0.4 kpc, B^ = 10 //G, and R" = 4- kpc. The 
strength of the halo magnetic field above the plane is shown in 
Fig.[II] 

We add the halo field to the disk field, run the simulation 
and obtain average longitude profiles for 25° < b < 35° and 
-35° <b < -25° as shown in Fig. [12] Some RMs of EGSs at 
longitudes of about 90° and larger show large deviations both 
above and below the plane. These areas coincide with the po- 
sitions of Loop II and Loop III. However, whether the Loop 
structures physically cause these RM deviations needs a further 
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Fig. 9. Different magnetic field configurations in the Galactic plane as seen from the north Galactic pole. The corresponding RM 
profiles along the Galactic plane (\b\ < 5°) are shown in the panel below. The average of RMs within longitude bins from the CGPS 
and the SGPS are shown for comparison. 
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Fig. 10. The same as Fig. |9l but for the ASS +RING, th e ASS+ ARM and the BSS field configur ations. Bes i de the binned CGPS 
and SGPS data the RMs near the center from lRov et al.l (l2005h and RMs from the compilation of iHan et alj ( Il999h are shown by 
triangles and squares, respectively. For the ASS+ARM model the NE2001 model is overlaid onto the field in gray scale. 
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Fig. 11. The strength of the regular halo magnetic field compo- 
nent. Contours start at 0. 1 fiG and increase by a factor of 2. The 
field below the plane has the same strength but is reversed. 



check. Except for these outliers the simulations are generally 
consistent with the observational data. 

We show the latitude variation of the simulated RMs in 
Fig. [T2l and in addition the average RM between / = 100° and 
/ - 120° versus latitude from -5° to 20°. The observed data 
are from the extended CGPS (Jo-Anne Brown et al., in prep). 
The halo field plus the ASS+RING or the ASS+ARM disk fields 
yield quite a good fit. However, the halo field plus the BSS field 
is clearly inconsistent with the data. If we tune the BSS field to 
match the latitude profile, we fail to get a good fit to the RMs 
along the Galactic plane. We conclude that the BSS model is not 
in agreement with the observed RMs of EGSs. 



6.2.3. Remarks 

We introduced field reversals in our models with sharp transi- 
tions, which means an abrupt field change in the opposite direc- 
tion. This might be not a realistic scenario, however, we do not 
expect much influence by this simplification as the transition re- 
gions are smoothed out in our simulated maps. We also used a 
constant pitch angle in our models, which is most likely another 
simphcifation. We note that the ASS+RING and the ASS+ARM 
models are not very sensitive to the pitch angle used. 



6.3. Random fields, CR electron density and local excess of 
the synchrotron emission 

6.3.1. The random magnetic field component 

The knowledge about th e properties of random magnetic fields in 
galaxies is quite limited ('Beck'2001'). A typical scale for the ran- 
dom fields wa s claimed t o be about 50 pc based on the single cell 
size model ( Rand & Kulkarnilll989h . However, the single cell 
model is a strong simplification. Structure functions of RMs of 
EGSs yield the fluctuation properties of rigbw. Unfortunately, due 
to the coupling between the electron density and magnetic field 
the extr acting of the magnetic fiel d fluctuation remains very un- 
certain. [MinteF&^^gangle^ ([199^ tried to tackle the decoupling 
by involving the structure function of emission measures and 
concluded that the magnetic field spectrum can be described by a 
Kolmogorov spectrum with an outer scale of about 4 pc. The tur- 
bulence for the region of sky investigated by Minter & Spangleg 
(1 19961) occurs in thin sheets or filaments and can be described as 
a two-dimensional turbulence for scales between 4 pc to about 
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80 pc. lHan et al.l (|2004|) claimed a flatter spectrum for scales be- 
tween 500 pc and 1 .5 kpc based on pulsar RM and DM data. 

In principle a stochastic realisation of the random magnetic 
fields can be generated by the Fourier transform of the square 
root of the spectrum as (e.g. lMartell l2005) b(s) cc J VP(k)(^ -H 
i^)exp(/k ■ s)ds, where P(k) is the power spectrum, k is the 
wavenumber, s is the spatial position, and ^ and ( are Gaussian 
random numbers. This requires a 3D inverse Fourier transfor- 
mation to obtain b(s). To fulfill the above equation the Galaxy 
must be gridded to small cubic boxes with a volume of P. To 
avoid an intrinsic correlation by the selected box size, they must 



be smaller than the outer scale of the magnetic field fluctua- 
tions. Let / have a size of 2 pc and the Galaxy volume to be 
40 kpcx40 kpcxlO kpc 2 x 10'^ boxes for an all-sky calcula- 
tion are needed, which is far beyond current computer memory 
capabilities. However, a small patch of the sky can be treated 
that way, utilizing a huge amount of computing time. A more 
detailed discussion on how to obtain a random field from a 
Kolmogorov spectrum for high-angular resolution simulations 
will be presented in a subsequent paper, where we will simu- 
late high-resolution patches of Galactic emission as they will be 
observed by LOFAR and the future SKA. 

Here we assume that the random field components are ho- 
mogeneous and follow a Gaussian distribution in strength with 
an average of zero and a scatter of cr^/ V3 in each dimension. As 
discussed later the best value for the mean random field strength 
is b = 3 juG throughout the Galaxy. 

6.3.2. CR electron density 

Synchrotron emission in the frequency range from 10 MHz to 
40 GHz originates from CR electrons with energies between 
400 MeV and 25 GeV assuming a magn etic field cornponen t 
perpendicular to the line of sight of 4 fiG (IWebber et al.lll980h . 
Local measurements of the density of CRs with energies below 
several GeV are very uncertain due to the limited know ledge of 
the heliospheric modulation (e.g. IStrong et alJl2004h . For the 
CR electrons w ith energies higher than several GeV, the mea- 
surements (e.g. iBarwicket al.lll998t iGrimani et aP 12002 !) fol- 
low a simple power law. The flux density can be described as 
N(y) = C(R,z)y'' with p of about -3, where y is the Lorentz 
factor. The CR flux density at 10 GeV from the aforemen- 
tioned measurements is about 0.2 (GeV m^ sr s)"' which means 
C(R = 8.5 kpc,z = 0) = Co ~ 3.2 x 10"^ cm^^ However, it 
is questionable wheather the local measuremen ts are representa- 
tive for the CR density el sewhere in the Galaxy (iPohl & Espositol 
119981 : IStrong et al.ll2004 ). Therefore we regard below the mea- 
sured local CR flux density as weakly constrained. 

We assume that the CR electron distribution consists of a 
disk component, whose spatial distribution is written as 



C(R,z)^Coexp 



R-Rp 



8 kpc 1 kpc 



(11) 



with Co = 6.4 X 10"^ cm"^. For the Galactic centre region the 
CR electron flux density is taken as a constant C(R < 3 kpc) = 
C(R = 3 kpc). The disk component is set to zero for z > 1 kpc 
as discussed in Sect. 7.2. 

The spectral index of the CR electrons p, the radio flux den- 
sity a and the brightness temperature /3 are related to each other 
as a = (p + l)/2 — p + 2. Spectral index maps show spatial 
variations of the spectra l index and a spec tral steepening at high 
frequencies is indicated jReich et al.l2004l) . In this paper we sim- 
plify the problem and keep the energy spectral index of the CR 
electrons p constant with a value of p - -1> corresponding to 
/? = -3. For the frequency range below 408 MHz we assume p 
tobe-2orj6 = -2.5. 

6.3.3. The local excess of the synchrotron emission 

The 408 MHz total intensity map shows a nearly constant high- 
latitude {\b\ > 30°) emission, which can be explained in two 
ways. There might exist an isotropic disk-centered halo emission 
component, which is difficult to model from our location in the 
Galactic plane far from the centre. However, the halo magnetic 
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Fig. 14. The same as Fig. [T3] but for profiles of polarized intensity at 22.8 GHz. All the maps are smoothed to an angular resolution 
of 2°. 



field we derived from RMs makes a smooth total intensity halo 
emission component quite unlikely. Otherwise the high-latitude 
emission can be attributed to local emission. In fact the observed 
local synchrotron enhancement requires a local excess of the CR 
electron density and/or an increase of the local random magnetic 
field component. 

Absorption measurements towards distant HII regions lo- 
cated in the inner Galaxy at 74 MHz yield a synchrotron emis- 
sivity of 0.35 K pc"' to 1 .0 K pc"' for distances larger than about 
2 kpc dNord e t al. 2006), which is consistent with the 408 MHz 
emissivity from the model of lBeuermann et akl (11985 ) for a typi- 
cal sp ectral index^ = - 2.5 for low frequency synchrotron emis- 
sion JRogeret al.lll999(). The 22 MHz survey absorption data 
from HII regions by 'Roge r etaLl d 19991) give an emissivity of 
about 30 K pc ' for distances less than about 1 kpc, which is 
a f actor of about 2 larger th an the typical emissivities obtained 
by iBeuermann et al.l (Il985l) . indicating a local synchrotron en- 



hancement by a factor of 2. [Fleishman & TokarevI d 1995b sum- 
marized low-frequency absorption measurements up to 2 kpc, 
which indicate a local synchrotron enhancement by a factor of 
about 3. From polariza tion modelling towards th e very local 
Taurus molecular cloud IWoUeben & ReichI (|2004|) confirmed a 
significant local synchrotron excess. 

The solar system is located in a low-density region filled with 
hot gas, called the "local bubble". Its size is about 200 pc. The 
"local bubble" is believed to be created by about twenty super- 
novae from the nearby Scorpius-Centaurus OB clusters during 
the last 10- 15 Myr. The last explosion may have occurred abou t 
1 Myr ago (i Breitschwerdt & de Avillezl2006l:lFuchs et alj2006h . 
There is compelling evidence that shock accel eration in SNRs 
is the source for energetic CR electrons (e.g. iKobavashi et al.l 
2004 ). With a t ypical diff'usion coefficient of lO'^* cm-^ s~^ 
(Ko bavashi et al.l 12004) CR electrons may have reached about 
600 pc distance and thus increase the CR electron density in the 
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local bubble. Alternatively, multiple fragments with compressed 
magnetic fields from the shells of the numerous very old SNRs 
forming the "local bubble" might have survived in the interstel- 
lar medium and enhance the local synchrotron emission. From 
this scenario the local excess of the synchrotron emission can be 
achieved either by an increase of the local CR electron density 
and/or by an enhanced local random magnetic field as akeady 
stated earlier We have tested both cases and find that the results 
are the same. For the simulation results shown below we assume 
a local spherical component within 1 kpc from the Sun, where 
the CR electron density exceeds the disk component by a factor 
1.5. 



6.3.4. Simulated 408 MHz total intensity and 22.8 GHz 
polarized intensity 

In Figs. [13] and [14] the latitude and longitude profiles from the 
simulated maps are compared with those from the 408 MHz and 
the WMAP 22.8 GHz polarized intensity all-sky surveys. All 
three models fit the total intensity profiles quite well. However, 
for polarized intensities the BSS magnetic disk field cannot well 
reproduce the small longitude asymmetry relative to the Galactic 
centre along the Galactic plane (Fig.[T4ll. 

The direction of the halo magnetic field is in the opposite 
direction to the disk field above and in the same direction below 
the Galactic plane. Therefore the resulting total magnetic field 
strength is stronger below than above the Galactic plane. In our 
model the CR electron density is truncated at z = 1 kpc and the 
halo field decreases quickly towards the Galactic plane, which 
makes the resulting intensity diff'erence quite small, e.g. about 
20 /iK for polarized 22.8 GHz emission (see Fig.[T4]i. Therefore 
this asymmetry is difficult to prove observationally and requires 
a careful absolute zero-level setting for the large-scale emission. 
For the case of a smaller halo magnetic field and an extent of the 
CR electrons beyond the present limit of 1 kpc in z this difference 
will remain small. 



6.3.5. Remarks on the random Galactic magnetic field 

We assume a random field component of uniform strength 
throughout the Galaxy. However, the random magnetic field 
could be considered to be stronger in the Galactic plane, where 
supernova explosions and violent star formation processes in- 
ject energy by shock waves or stellar winds into the interstel- 
lar medium and thus initiate higher fluctuations of the magnetic 
field than at high latitudes. From the analysis of the structure 
function of EGSs it was shown that fluctuations of the magne- 
tized inters tellar medium in t he Galactic plane are longitude- 
dependent JSun & Hani |2004|) and also d iff'erent for arm and 
inter-arm regions (JHaverkorn et al.l [20061) . Our actual simula- 
tions are not sensitive to spatial variations of the random mag- 
netic field at high-latitudes because of the present truncation of 
the CR electrons at z = 1 kpc. 



6.4. Depolarization 

Using the Galactic magnetic field model and the CR electron 
distribution we are able to model the observed total intensity 
map at 408 MHz and the polarized intensity map at 22.8 GHz. 
The 1 .4 GHz total intensity map has a different amount of ther- 
mal and non-thermal components and is also well modelled (not 
shown here). The 1.4 GHz polarized intensity originates from 
synchrotron emission modulated by the interstellar medium and 
should be reproduced as well after accounting for depolarization 
effects, which needs to be discussed in some detail. 



6.4.1. Simulation problems 

We show the simulated results for polarized intensity at 1 .4 GHz 
in Fig. [15] where the longitude profiles along the Galactic plane 
and the latitude profiles for |/| < 30° are presented. The ob- 
served profiles show strong depolarization observed in partic- 
ular towards the inner Galaxy (see Fig. |2]i. Note that promi- 
nent large structures like the "Fan region" and Loop I (as in- 
dicated in Fig.[T5]l are not included in our modelling, as they are 
considered as local structures. Therefore the lower envelopes of 
the observed profiles have to be compared with the simulations. 
Towards longitudes / - ±90° the magnetic field direction is al- 
most tangential to the line of sight, which makes any polarized 
intensity intrinsically small and the simulated polarized inten- 
sities at 1.4 GHz are almost consistent with the observed ones. 
However, towards the centre and the anti-centre region the sim- 
ulated polarized intensity is in general much larger than that ob- 
served. 

We show in the following that the random magnetic field 
component is properly set in our simulation: Assuming that the 
intrinsic random field component is located in cells with a com- 
mon size /, the number of cells in the i-\h volume unit calculates 
then as A^^ = fiAr^A6^/P. In the outer (large /) unit the number 
of cells is very large. However, in the simulation we specify also 
for the outer unit just one direction for the random magnetic field 
component, which means the number of random magnetic field 
cells is clearly underestimated. This way the simulated polarized 
emission may be too large. If this is the case the simulated polar- 
ized emission should significantly decrease when the maps are 
smoothed. We smoothed the simulated polarization map to an 
angular resolution of 36', which is the resolution of the observed 
1.4 GHz all-sky map, and to a resolution of 2°. The results are 
shown in Fig.[T6]for a cut along the Galactic plane. It is evident 
that the changes in the profile after smoothing are small, which 
indicates that the cell size for the random field component is 
small enough even for large radii. This proves that the random 
field component to be properly considered in the simulations. 



6.4.2. Solutions to increase the depolarization 



The local total magnetic field vector Btot = B-i-b in our model 
is about 4 yuG with a ratio between the random to the regular 
magnetic field component of about 1.5. We note that this total 
field strength is just consistent with the value of 6+2 yuG, which 
was derived by assuming energy equipartition between CRs and 
magne tic fields from the synchrotron model of Beuermann et al. 
( Il985l) (Berkh uijsen, personal communication, shown in Fig. 1 
in iBeckllMTl) . 



As shown by Sok oloff et al.l (1 1 9981) Faraday depolarization ef- 
fects become strong when either the RM or the RM fluctuations 
along the line of sight (ctrm) increases. For our simulations the 
RM fluctuations can be estimated: a-f^M ~ 0.81 < «(. > cri, VrAr, 
where < n^ > is the average thermal electron density. In our 
calculation Ar = 20 pc. For typical values such as r = 10 kpc, 
< rie >= 0.03 cm"-', crb - 3 yuG, the RM fluctuation can be 
estimated to be about 32 rad m"^. This fluctuation should pro- 
duce enough depolarization at 1.4 GHz based on the formula 
by Sokoloff et al. (1998). However, the assumption of identical 
distributions of all constituen ts of the magnetic-io nic medium 
along the line of sight made bv lSokoloff et alj (Il998l) clearly dif- 
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Fig. 15. Simulated polarized intensities at 1 .4 GHz (shown as 
thick solid and dashed lines) along Galactic longitude in the 
plane (upper panel) and along Galactic latitude (lower panel) are 
shown in comparison to the 1.4 GHz observations (thin solid 
lines). The simulation shown in dashed lines includes a correc- 
tion for RM fluctuations based on the thermal electron filling 
factor and a related magnetic field enhancement. For details see 
text. The angular resolution is always 36'. 



where flie factor f„i, is defined as 



fnb = 



< tie >^ (tI 
< nlb^ > 



(13) 



Then we calculate the RMs in each volume unit as RM - 
RM() + RMrl sjf^, where RMq and RM,- are calculated from 
the regular and random magnetic field components, respectively. 
From the simulations we found that for f„h - 0.01 the polarized 
intensity in the Galactic plane is significantly reduced and the 
simulated results are consistent with the observations, as seen in 
Fig. [15] The correction does not change the total RM along the 
line of sight on average, which we prove in Fig. [17] because just 
the random RM component in each volume unit is increased. 



K 




180 



1(°) 



Fig. 17. The original simulated RM distribution for the 
ASS-hRING disk field (thin line) and the RM distribution in- 
cluding RM fluctuations by fnb = 0.01 (solid line). 




180 



Fig. 16. Simulated polarized intensities at 1 .4 GHz along the 
Galactic plane are shown at 15' (thin solid line) and after 
smoothing to 36' (thick solid line) and 2° (dashed line). 




Fig. 18. The spatial fluctuations of RMs before (dashed) and 
after (solid line) considering RM variations (/„/, = 0.01) in a 
clumpy magnetized WIM along the line of sight. 



fers from the more realistic properties used for our simulation. 
The observed depolarization at 1.4 GHz suggests that this RM 
fluctuation is too small to cause significant depolarization. 

The mean RM has already been constrained by extragalactic 
sources. We try to increase the RM fluctuations in each volume 
unit to introduce a larger depolarization. We rewrite the formula 
for RM fluctuations for the case of a clumpy thermal medium 
using the same units as in Eq. (1), which reads; 



RM 



0.81^ < nlb^ > rAr -■ 



0.81' 



, < n, 



>'-l 



fnb 



rAr, 



(12) 



What is the physical meaning of the factor f„h ? If the ther- 
mal electron density and the random field are independent, i.e., 
< nlb^ >-< nl >< b^ >, the factor f„b is identical to the ther- 
mal electron filling factor fg. Note that < b^ >- cr^. The ther- 
mal electron fillin g factor in the pla n e is a bout 0.07 from our 
simulations, while iBerkhuiisen et al.l (120061) derived values for 
/(. between 0.04 and 0.07, which is larger than the value /„;, = 
0.01 as obtained above. This indicates a coupling between the 
thermal electron density and the random magnetic field strength. 
In general we may write f„h - fefc, where we take fc as the 
coupling factor In case a dependence like b^ oc («^)" exists. 
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we have fc oz f^. Our simulations clearly favor /,;, < /, with 
a > Q. This means that the thermal electron density and the ran- 
dom magnetic field strength are correlated. However, we note 
that the physics and a theoretical value for the coupling factor 
needs to be worked out. Therefore we take our results as qualita- 
tive and did not attempt to account for possible spatial variations 
of /„/, as it is suggested by /^(z). 



6.4.3. Comment on the observed RM fluctuations 

We finally discuss the observed large fluctuating RMs from 
EGSs in the Galactic plane (see Fig. |4]i. Note that these fluctu- 
ations are in principle different from those we discussed above. 
RM variations of EGSs result from different ISM properties for 
each line of sight, while RM fluctuations causing depolarization 
are along a single line of sight. However, larger RM variations 
are expected for different line of sights, if the RM fluctuations 
are large at all. Therefore both cases are related. As shown in 
Fig.[T8]the spatial RM fluctuations are significantly larger after 
including the factor /„& = 0.01. We note that the large fluctu- 
ations are now consistent with the RM observations presented 
in Fig. |4] An important implication is that the large spatial RM 
scattering is not due to fluctuations of large random field com- 
ponents, but to the thermal filling factor and a coupling factor 
between the electron density and the random field strength of 
the clumpy magnetized WIM. 



7.2. The truncation of the CR disk electrons at 1 kpc 

CR electrons gyrate in a magnetic field with a radius rg =s 
1 X lO^^E/B pc, where the energ y E is in units of GeV and the 
magnetic field strength is in fiG (lFerrierell200ll) . Turbulence in 
the halo can be inferred from super-bubbles generating shock- 
waves from supernovae in OB clusters and stellar winds. The 
outer scale can be estimated to be about 100 pc, which is much 
larger than the gyrating radius rg. In that case diffusion of CRs 
both along and across magnetic field lines is strong. The dif- 
fusion along magnetic field lines is much larger when the tur- 
bulent component is smaller than the regular field component 
( Casse et al. 2001) as it is the case for the present halo model. 
Although the d iffusion of CRs in a turbulent field is still debated 
(|Ptuskinll2006) . it seems unlikely that the CR electrons can be 
prevented to diffuse out of the Galactic plane into the region with 
the large halo field. Therefore the truncation of the CR electron 
density at 1 kpc we need for our simulations seems unrealistic. 

However, if we give up the 1 kpc CR electron truncation the 
modelled synchrotron emission will increase largely at high lat- 
itudes. As a consequence two ridges of strong emission running 
parallel to the Galactic plane will show up, which are not ob- 
served. In case we reduce the halo field to avoid the synchrotron 
emission problem, we are unable to model the RM data at high 
latitudes. To avoid this it is required to increase the thermal elec- 
tron density or the scale height of the DIG. 



6.5. Summary 

We summarize the modelling constrains and the related obser- 
vations in Table [1] and show the simulated 408 MHz total inten- 
sity, the 1.4 GHz polarized intensity, 22.8 GHz polarized inten- 
sity and polarization angle at 22.8 GHz in Fig.[T9]based on the 
ASSh-RING disk field configuration. All maps shown have an 
angular resolution of 15' and can be obtained on request from 
the authors. 



7. Discussion 

7.1. Is a large halo field physically reasonable? 

Halo magnetic fields were detected just in a handful of nearby 
galaxies such as NGC 891 with a halo field running parallel 
to the plane or NGC 4631 showing a vertical field direction. 
The field strength i s between 4 uG an d 6.5 fxG for z between 
2 kpc to 5 kpc (e.g. lDumke et alj|1995h . However, the observed 
high-latitude RMs require a stronger halo field for the Galaxy as 
shown above. 

In case the dynamo theory is relevant to describe the Galactic 
magnetic field evolution we learn that in a thick disk a asymmet- 
ric hal o field can be ge nerated as shown by Soko loff & Shukurov 
( Il990l) . Simulations bv lBrandenburg et alj ( IT992 ) indicate that it 
takes a long time, which is comparable to the Hubble time, that 
the halo field reaches stability. 

Considering hydrostatic equilibrium in the halo a mid-plane 
magnetic field strength of about 4 fiG with a scale height of 
4.4 k pc in the solar neighbourhood is needed JKalberla & KerpI 
11998 ). The present halo model gives a field strength of about 
7 fiG at z of 1.5 kpc. This field sharply decreases towards the 
Galactic plane. Towards larger z our model is insensitive to any 
scale height. 



7.3. The scale height of the DIG 

The scale height of the DIG of about 1 kpc results from fitting 
DM sin 1^1 versus \z\ for pulsars with in dependent distance mea- 
surements (iGomez et al.l 120011: ICordes & Lazio 2002) . In fact 
the most stringent constraints come from distant pulsars located 
in globular clusters. However, there are just very few clusters 
with a large z and the ISM has many voids and clumps, which 
may influence the scale height derived from limited data. 

From the WHAM Hg survey observations in the Perseus arm 
direction dHaffner et al.lll999l: [Reynolds etakll 19991) the Ho- in- 
tensity can be fitted with a scale height of 1 kpc and a mid-plane 
density of 0.06 cm""*. This electron density is larger than the 
mid-plane density of about 0.03 cm"^ in the NE2001 model. 
Since the NE2001 model fits the DM sin |/7|-|z| relation very 
well, this might indicate that the electron density at high z could 
be larger The WHAM survey clearly shows WIM emission 
enha ncements for large |^| relative to the EM sm\b\ distribu- 
tion JHaflrner et al.lll999HRevnolds et alJI 19991) . which indicates 
the presence of areas of large electron densities at high lati- 
tudes. However, how these high-latitude features are distributed 
is rather unclear at present. 

In order to check how the high-latitude electron distribution 
influences the present simulations we run the HAMMURABI 
code by increasing the scale height of the DIG to 2 kpc by keep- 
ing the mid-plane electron density from the NE2001 model un- 
changed. As a consequence from this model a regular halo field 
magnetic strength of about 2 fiG with a scale height of 4 kpc is 
required, which sounds more reasonable compared to the 10 fiG 
magnetic field strength enforced by a 1 kpc DIG scale height. 
In this case the halo field param eters agree with the estimates 
of ISokoloff & Shukurovl (1 19901) based on the dynamo theory. 
The 1 kpc truncation of CR electron distribution could also be 
avoided that way. We conclude that the high-latitude thermal 
electron density is a critical issue for Galactic 3D-modelling and 
requires more investigations. 



18 X.H. Sun et al.: Galactic 3D-models 

Table 1. Observational features versus model constraints based on the NE2001 thermal electron model. 



observational features 



model constrains 



22.8 GHz free-free emission template andj822/408 mhz 

EGS RM profile along the Galactic plane 

High-latitude RMs asymmetry (to the plane and center) 

RM latitude profile at 100° < / < 120° 

Longitude PI asymmetry to the centre at 22.8 GHz 

Synchrotron emission longitude profiles: 408 MHz 1/22.8 GHz PI 

Isotropic high-latitude emission at 408 MHz 

Strong depolarization at 1 .4 GHz 



NE2001 plus filling factor for thermal electrons 

ASS+RING, ASS-hARM or BSS regular disk field, local field strength 2 yuG 

Asymmetric toroidal halo field, field strength up to 10 yuG, CRs up to 1 kpc 

Favour ASS+RING or ASS-hARM model 

Favour ASS+RING or ASS+ARM model 

Disk CR electron component and 3 fiG random field component 

Local excess of the synchrotron emission 

Large RM fluctuations along the line of sight 

by a small n^ filling factor and a coupling of the random field with n^ 



I 4^^m»£ FSSrlUNG 
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Fig. 19. Simulated all-sky maps for total intensity at 408 MHz, polarized intensity at 1.4 GHz and at 22.8 GHz and polarization 
angle at 22.8 GHz are shown for the ASS+RING magnetic field configuration. All the maps have a resolution of 15'. 



8. Conclusions 

Galactic 3D-models of the distribution of thermal electrons, 
cosmic -ray electrons and magnetic fields were used to calculate 
via the HAMMURABI code all-sky total intensity, polarization 
and RM maps. A comparison with observed maps over a wide 
frequency range is made to adapt the 3D-model parameter in the 
best possible, although still qualitative way (Table 1). 
The main conclusions are : 

- The NE2001 thermal electron model combined with a filling 
factor /e as derived by Berkhuiisen et al. (2006) reproduces 
the optical thin 22.8 GHz free-free emission template and 
also the observed spectral flattening in the Galactic plane be- 



tween 408 MHz and 22 MHz caused by absorption of op- 
tically thick thermal gas. We find indications that the fill- 
ing factor might be large r towards high z than that found by 
iBerkhuiisen et al.1 (12006'). 

An ASS field configuration plus a reversal inside the so- 
lar circle either in a ring or between the inner edges of the 
Sagittarius-Carina arm and the Scutum-Crux arm is favoured 
by RMs of EGSs observed along the Galactic plane. A BSS 
field configuration also fits these RMs, however, it then fails 
to fit the observed RM gradient in latitude direction. The 
strength of the regular magnetic field component at the so- 
lar radius is 2 fiG. The field in the Galactic centre region is 
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constant at 2 jjG in a direction opposite to that at the solar 
radius. 

- The strength of the Galactic random field component is 3 fiG. 

- A regular halo magnetic field with a maximum of 10 fiG at 
a height of 1.5 kpc is required to account for the observed 
high-latitude RMs. This toroidal field is asymmetric to the 
Galactic plane and decreases sharply towards the plane. The 
variation of the field towards high z is not a sensitive param- 
eter in our simulation. 

- The CR electron distribution has an exponential scale height 
of 1 kpc and a radial scale length of 8 kpc. The disk is 
truncated at a height of 1 kpc. The local CR flux density at 
10 GeV is taken twice the directly measured value at earth. 
The energy power law spectral index p is taken to be constant 
at -3 for frequencies above 408 MHz and -2 below. 

- The observed local excess of synchrotron emission is taken 
into account. The excess is either caused by an enhanced lo- 
cal CR electron density or by an increase of the local tur- 
bulent magnetic field. Our modelling can not distinguish be- 
tween both possibilities. 

- In addition to the filling factor we introduce a coupling factor 
between the thermal electron density and the random mag- 
netic field component to reproduce the observed large RM 
fluctuations and the observed depolarization at 1 .4 GHz. 

However, a number of caveats remain to the above conclu- 
sions. In our simulation we use the NE2001 model as input 
for the thermal electron distribution and density. This requires 
a large halo field strength and hence a truncation of the CR elec- 
trons at z = 1 kpc. This might be not very realistic. Indeed, there 
are indications as discussed in Sect. 7.3 that the scale height of 
thermal electrons might be larger than 1 kpc. When the thermal 
electron scale height is increased by a factor of 2 we get a regular 
halo magnetic field of 2 juG and the truncation of the CR electron 
density at 1 kpc is avoided. We conclude that the high-latitude 
thermal electron density needs more investigations and thus the 
halo field parameters remain open. 

For large areas along the Galactic plane RM data are still 
missing or very rare. A final answer on the magnetic field config- 
uration depends on a much more complete RM data set. Pulsar 
RM data are not yet included in our model, although they are 
potentially better tracers of the regular Galactic magnetic field 
along the line of sight than those from EGSs as they show no 
intrinsic RM. Pulsar RMs are clearly important to detail the disk 
magnetic field throughout the Galaxy on smaller scales than in 
our model and might also reveal deviations from a general reg- 
ular pattern being assumed here. Despit e of significant obser- 
vational progress over the last years (e.g. Han et all2006 ^ more 
RM data are clearly needed to be conclusive and pulsar distances 
need to be better constrained in addition. 

We note that our 3D-model diff'ers from previous models, 
which were derived from selected observational data. These 
models then fail to properly reproduce other observations, 
which, however, trace relevant information of the interstellar 
medium, which can not be ignored on the way towards a physical 
relevant model. 
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